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PROJECT  SUMMARY 


The  purpose  of  this  research  program  was  to  expand  our  understanding  of  the 
growth  of  fibers  of  the  nonlinear  optical  material,  MgO.UNbO<.  This  material  is 


or  interest  because  of  its  large  nonlinear  optical  coefficients  and  its  good 
optical  properties  The  undoped  material  suffers  from  photorefractive  damage 
It  nas  been  found  that  by  doping  UNb()3  with  approximately  5%  of  Mgn  the 


Dhasematch  temperature  increases  to  the  point  that  the  material  becomes 


self  annealing.  This  makes  this  material  particulaly  attractive  for  second 
harmonic  generation  of  infrared  laser  sources  and  for  other  nonlinear  devices. 
Second  harmonic  generation  of  laser  sources  usinq  fiber  materials  has 


several  important  advantages.  Such  components  could  be  used  In  integrated 
optical  devices,  significantly  reducing  the  size  of  these  systems.  Moreover,  the 
use  of  fibers  for  these  elements  can  lead  to  Important  improvements  in  the 
efficiency  of  conversion  of  the  laser  source  to  the  second  harmonic  as  the  guided 
wave  will  maintain  collinearity  over  the  length  of  the  fiber.  This  increases  the 
interaction  length,  resulting  in  an  improved  conversion  efficiency.  Fibers  are 
also  very  easy  to  temperature  control,  an  important  feature  for  phase  matched 
second  harmonic  generation. 

During  the  course  of  our  project  we  studied  the  growth  characteristics  of 
MgO.LiNbOx.  We  made  18  growth  runs  of  this  material  and  were  able  to  improve 


nn  the  growth  conditions  so  as  to  be  able  to  grow  fibers  of  good  optical  quality 
without  resorting  to  post  growth  processing  However,  the  mechanical  strength 
of  toe  fibers  were  not  as  desired  and  an  annealing  process  was  required  to 
increase  the  strength  of  the  as-grown  fibers.  Problems  were  encountered  when 
we  attempted  to  pel?  the  two  domain  a-axis  fibers  in  order  to  obtain  single 
domain  fibers  Asymmetric  heating  and  electric  field  poling  and  both  techniques 
together  were  tried  with  out  success  Control  ofthe  process  was  too  critical  for 
the  system  we  used  ann  several  fibers  were  lost  The  only  process  chat  appears 
premising  to  us  is  symmetric  heating  with  electric  field  polmg  near  the  melting 
pmnt  ot  figi  >  ijNb(  w  Further  work  needs  to  be  done  to  evaluate  the 

characteristics  of  this  process  in  order  to  obtain  single  domain  fibers  of  this 
material  Once  this  is  achieved,  highly  efficient  nonlinear  optical  devices  using 
fibers  of  MqO.LiNbO^  will  result 


SECTION  1 


INTRODUCTION 


As  spectroscopic  and  other  optical  instruments  are  used  in  a  wider  range 
of  experiments  and  environments,  such  as  in  outer  space,  it  becomes 
necessarv  to  make  these  instruments  smaller  and  more  compact,  lighter  in 
weight  and  more  efficient.  The  most  efficient  and  compact  laser  source  at 
this  time  is  the  semiconductor  laser.  This  laser  has  achieved  efficiencies  of 
greater  than  50%.  However  this  laser  has  some  serious  limitations  that  do 
not  allow  its  use  for  many  scientific  studies.  Some  of  these  limitations  are 
that  the  beam  quality  of  these  lasers  is  usually  poor,  the  beam  divergence  is 
Great  and  the  beam  is  astigmatic  and  the  laser  can  not  be  tuned  over  a  wid^ 
spectral  range.  However  this  laser  is  ideal  as  a  pump  source  for  a  solid  state 
laser  that  can  overcome  these  shortcomings.  A  system  that  could  be 
envisioned  as  a  spectroscopic  source  is  a  semiconductor  diode  laser  pumping 
a  Nd.YAG  laser  that  is  frequency  doubled  to  532  nm.  This  output  is  then  used 
to  pump  a  T1.AI2O3  laser  that  is  continuously  tunable  from  660  nm  to  over 

1000  nm  The  possibility  thus  exists  for  an  all  fiber  laser  system.  Fiber 
lasers  of  Md  VAG  have  already  operated  successfully  The  next  step  would  be 
to  develop  the  frequency  doubler  to  double  the  output  of  this  fiber  laser  to 
532  nm  Significant  success  has  already  been  achieved  at  Stanford  University 
in  the  growth  of  single  crystal  fibers  of  MgOUNbOj  as  frequency  doublers. 

Laser&emcs  has  the  exclusive  license  to  the  patents  held  by  Stanford 
University  on  the  laser-heated  pedestal-growth  (LHPG)  technique  for  growing 
single  crystal  fibers.  We  have  investigated  the  growth  of  MgO.LINbCh  using 
this  technique  as  a  function  of  the  melt  temperature  and  the  atmosphere 
using  source  rods  prepared  from  single  crystal  material  grown  by  Crystal 
Technology,  Inc  We  also  investigated  the  post  growth  annual  to  minimize 
defects  and  oxygen  loss  in  the  fibers  We  have  made  preliminary 
investigations  ot  techniques  to  pole  these  fibers  in  order  to  obtain  single 
domain  a-avis  fibers  with  limited  success 

Mgr'i  l.)Mbf>3  fibers  are  potentially  useful  in  a  variety  ot  electo-opticai 

devices  that  take  advantage  of  its  large  piezoelectric,  pyroelectric, 
electro-optic  and  nonlinear  optical  coefficients.  Nonlinear  interactions  in 
the  fiber  form  are  attractive  because  the  fiber  eliminates  beam  spreading  due 
to  diffraction,  thus  leading  to  enhanced  frequency  conversion  efficiencies. 

Our  primary  interest  is  in  frequency  doubling  the  1.06  micron  radiation  from  a 
fiber  NdYAG  laser.  At  this  frequency  lithium  niobate  will  be  non- critically 
phasematched  near  room  temperature. '  To  obtain  non- critical  phasematching 
m  the  fiber,  the  fiber  must  be  grown  orthogonal  to  the  c-axis  It  has  beer, 
found  that  MqO.LINbQz  will  grow  readily  along  the  a-axis  Two  complications 


result  during  the  growth  of  MqO.LiNbO-<  fibers.  These  are  the  oxygen  loss 
during  growth^  and  fhp  poor  mechanical  ?trAnqth  of  the  fibers  Associated 
with  t.hp  oxygen  loss  Is  a  brownish  coloring  that  the  fibers  evhiDit,  making 
them  unsuitable  for  optical  applications  However,  it  is  possible  that  the 
loss  of  oxygen  can  he  corrected  by  annealing  the  fibers  in  a  wet  oxygen 
atmosphere  for  eight  hours  at  900°  C .2*3  This  annealing  procedure  leads  to  an 
additional  benefit  in  that  a  marked  improvement  in  the  mechanical  properties 
of  the  fibers  results  as  well,  presumably  due  to  the  removal  of  crystal 
strains  in  the  fibers 

Another  problem  to  be  overcome  in  the  development  of  these  fibers  as 
frequency  doublers  is  to  pole  the  fibers  into  single  domains.  As  is  typical  of 
ferroelectric  materials,  as-grown  crystals  tend  to  break  up  into  a  number  of 
regions  of  reversed  spontaneous  polarization  (ferroelectric  domains).  Since 
the  sign  of  most  of  the  material  parameters  differs  in  antiparallel  domains, 
the  efficiency  of  many  devices  depends  on  the  preparation  of  "poled"  fibers, 
i.e  fibers  containing  only  a  single  ferroelectric  domain. 

As-grown  bulk,  crystals  of  MgOiiNbOj  tend  to  exhibit  multi-domain 


structures.  These  domain  distributions  are  conventionally  explained  in  terms 
of  a  balance  between  depolarization  energy,  which  favors  a  highly 
multidomain  configuration,  and  domain  wall  energy,  which  favors  a  single 
domain  crystal.4  In  c-axis  rods  these  domains  tend  to  take  the  form  of 
concentric  rings,  while  in  a-axis  rods  the  domains  form  a  series  of  layers 
perpendicular  to  the  rod  axis.  The  domain  wall  spacing  that  minimizes  the 
sum  of  the  depolarization  energy  and  the  domain  wall  energy  in  an  a-axis  rod 
can  be  calculated.  The  result  for  the  domain  wall  spacing,  0,  as  a  function  of 
the  rod  radius,  p,  is  given  by; 


Q  =  vrp[(TT3/2  IO)e0(e  +  i )  a/Ps2]'/2 


where  e  is  tne  dielectric  constant  of  the  rod,  a  is  the  domain  wall  energy  and 
P5  is  the  spontaneous  polarization  per  unit  volume,  Because  the  material 

properties  are  not  well  known  near  the  Curie  temperature  and  since  the 
quantity  in  the  square  brackets  is  independent  of  the  rod  diameter  it  is  best 
to  estimate  the  domain  wall  spacing  in  the  fibers  by  scaling  the 
measurements  in  the  bulk  crystals  by  the  square  root  of  the  diameter.  The 
typical  domain  size  for  centimeter  diameter  bulk  a-axis  boules  is  several 
hundred  microns.  Therefore  we  would  expect  that  the  domain  spacing  in  a  500 
um  fiber  would  be  of  the  order  of  100pm. 

However,  it  has  been  observed  that  the  fibers  do  not  break  up  into  small 
domains  In  fact  the  c-  axis  grown  material  is  single  domain  and  the  a-  axis 
fibers  contain  two  domains  that  are  anti-parallel  with  the  wall  perpendicular 
to  the  c  axis.  These  two  cases  are  depicted  in  figure  1.1.  It  is  clear  then 


that  for  the  fibers  depolarization  and  domain  wall  energies  are  not  the 
complete  picture.  Fejer-  has  determined  that  the  thermoelectric  effect  can 
produce  electric  fields  large  enough  to  pole  the  fibers  in  this  way  during  the 
growth  process.  This  electric  field  produced  by  the  thermoelectric  effect  is 
given  by 


E  =  Q  VT 


where  E  is  the  electric  field,  T  is  the  temperature  and  Q  is  the  second  rank 
thermoelectric  power  tensor.  Near  the  Curie  temperature  Q  has  the  value  of 
0.8  mWK  in  MgOliNbO^  Because  of  the  huge  temperature  gradients  that 

exist  during  tne  finer  growth  process,  large  enough  fields  exist  to  pole  tne 
fibers  as  observed  For  the  c-axis  growth  case  the  major  temperature 
gradient  is  along  the  c-axis  so  that  tne  fiber  is  poled  along  this  axis  For  an 
a-axis  growth  finer  tne  major  temperature  gradient  is  radial  and  snii  along 
tne  c-axis  This  radial  dependence  along  tne  c-axis  leads  to  tne  nead-to-neao 
domains  observed  in  a-axis  growth  fibers 

Tnese  observations  suggests  a  technique  to  pole  the  a-axis  fibers  in  order 
to  oDtain  a  single  domain  finer  By  asymmetrically  heating  tne  finer  on  tne 
surface  where  tne  c-axis  emerges  from  tne  fiber  Put  not  on  tne  opposite  face, 
tms  temperature  gradient  can  cause  tne  fiber  to  spontaneously  pole  Another 
procedure  would  be  the  application  of  an  electric  field  along  the  c-axis  The 
most  effective  poling  procedure  may  be  a  combination  of  these  techniques 
we  nave  investigated  tne  growth  of  MgOtiNbOj  in  finer  form  using  the 

laser-heated  pedestal-growth  technique.  In  section  2  we  discuss  the  theory 
of  nonlinear  effects  in  fibers.  In  section  3  we  discuss  the  laser-heated 
pedest3l-growth  technique  and  why  it  is  so  attractive  for  growing  high 
quality  single  crystal  fibers  and  in  section  4  we  discuss  the  results 


SECTION  2 


THEORY  OF  NONLINEAR  FIBER  DEVICES 

In  this  section  we  will  consider  those  effects  that  play  an  important  part  in 
determining  the  efficiency  for  second  harmonic  generation  (SHG)  in  fibers  We 
begin  by  relating  the  well-known  results  of  SHG  in  bulk  materials.  This  was 
first  discussed  by  Armstrong,  Bloembergen,  Ducuing  and  Pershan-  in  1962 
starting  from  Maxwell's  equations.  These  equations  can  be  solved  In  closed 
form  for  the  plane  wave  case.  A  nonlinear  polarization  term  is  introduced  at 
each  of  the  frequencies  of  the  waves.  These  nonlinear  polarization  terms  are  of 
the  form 


PmL\id- )  -  dt E(o)j )  Efcty )  2.1 

For  second  harmonic  generation:  coj  =  co |  +  to^  and  coj  =  coy  so  that  co(  =  2coj 

if  it  is  assumed  that  there  are  no  absorption  or  other  losses,  that  neglibie 
power  is  transferred  between  waves  and  that  the  amplitudes  vary  slowiy,  the 
wave  equation  can  be  easily  integrated  to  give  the  beam  intensity  at  the  second 
harmonic  as. 

1(2©')  =  to2  K2  deff2  L2  i2(0)  smc2(AkL/2)  2  2 

where  der f  is  the  effective  coupling  coefficient  which  is  related  to  the 
nonlinear  susceptibility  ana  Ak  is  given  by- 

At  =  2ka>  -  k.2w  =  2to  4 (e0M0)[n(«)  -  n(2©)]  2  3 


From  this  we  can  see  that  the  conversion  efficiency  increases  both  with  the 
interaction  length  L  and  with  the  intensity  1(0)  of  the  pump  beam.  One  would 
think  that  the  conversion  efficiency  could  be  increased  bv  focusing  the  pump 
beam  to  obtain  a  higher  value  of  1(0).  However,  the  effects  of  focusing  on 
conversion  efficiency  can  be  determined  bv  recalling  that  a  Gaussian  beam  of 
wavelength  a  focused  to  a  radius  w0  in  a  solid  with  refractive  index  n  remains 
collimated  over  an  approximate  length  2Cp,  where  the  Rayleigh  length,  Cr,  is 
given  by 

tp  =  Trnw0/'X  2.4 

Therefore,  as  the  beam  is  focused  with  shorter  and  shorter  focal  lengths  to 
decrease  w,v  the  interaction  length,  2tp,  is  decreased  by  the  square  of  wp.  With 


a  nonlinear  r  ceraction  involving  focused  beams,  a  trade-off  exists  between  a 
small  focai  spot  and  a  long  interaction  length.  For  noncritically  phase-matched 
second  harmonic  generation  in  a  crystal  of  length  L,  it  can  be  shown  that 
optimum  conversion  occurs  for  confocal  focusing,  i.e  ,  when  the  confocal 
parameter  of  the  pump  beam  is  approximately  equal  to  the  length  of  the  crystal 
Mathematically,  if  the  waist  of  the  Gaussian  beam  is  w0,  then  the  optimum 

ronjcino  condition  is  when 


L  -  5  68  2.5 

When  this  focusing  condition  is  met,  ne- conversion  efficiency  for  bulk 
second  harmonic  generation  in  the  undepleted  pump  limit,  is; 

=  I  (2o>)/ 1  (to)  =  (2co)2/e0p03  (d2/2n3)l(o)(  1  2L/X)  2.6 

Trie  conversion  efficiency  expression  for  guided  waves  has  also  been 
calculated  to  be  given  by 

tjc  =  K2(okM(a))  =  ( 2<o > ^ |jl q 3  (d2/2n~>)[l(co)/Aef f ]L2  L-‘ 

If  we  take  the  ratio  of  the  guided  and  bulk  efficiencies  from  Eqs.  2.G  and  2.7 
we  note  that 

gfv  T|p,  -  LX/(2.4nAeff)  *xG 

!t  is  interesting  to  determine  the  significance  of  this  result  We  define  a 
dimensionless  effective  area,  A.  by: 

A  =  Aeff/(TTpz)  2.9 

where  p  is  the  fiber  radius  Then  the  ratio  of  the  guided  to  bulk  efficiency  is 
Given  Pv 

J  t 

=  >LX/p^)(2.4nnA)"  *  2  10 

From  this  we  can  see  that,  for  a  given  radius,  the  relative  efficiency  of  the 
guided  wave  interaction  increases  with  both  L  and  X,  as  one  would  expect  as  the 
important  of  diffract  ion  alsn  increases  with  both  L  and  X  The  improvement  in 
efficiency  can  be  substantial  If  we  let  L=!0cm,  X  =  lpm,  p=25pm,  n=2,  A=l/2, 


6 


then  tiryiift  ~  40  We  have  developed  this  for  second  harmonic  generation  but  one 


would  evpect  similar  improvements  for  other  nonlinear  interactions  For 
interactions  involving  widely  separated  frequencies,  even  larger  improvements 
are  to  he  expected  because  of  the  relatively  poor  overlap  for  bulk  interactions 

For  a  lOcm-long  piece  of  lithium  mob3te,  the  efficiency  for  bulk  confocal 
second  harmonic  generation  at  532nm  is  A%! W  In  a  fiber  with  the  same 
parameters,  the  conversion  efficiency  would  be  .  16%/mW.  At  this  level,  one 

could  efficiently  double  weak  CW  sources,  eg.,  4mW  of  green  light  could  be 
obtained  from  a  50  mW  infrared  source.  The  gain  for  a  parametric  oscillator 
operated  at  degeneracy  is  equal  to  the  harmonic  conversion  efficiency,  thus 
milliwatt  thresholds  for  such  devices  would  seem  possible  as  well. 

In  our  analysis  to  this  point,  we  have  assumed  that  the  interaction  is 
noncntically  phasematched,  and  that  the  fiber  interaction  is  between 
fundamental  modes  with  the  crystal  axes  aligned  with  the  fiber  axes.  In 
critically  phasematched  bulk  interactions,  the  Poyntinq  vector  walkoff  limits 
the  minimum  useful  spot  size  for  a  given  length  of  crystal.  This  severely  limits 
the  conversion  efficiency  and  decreases  the  quality  of  the  output  beam  in 
fibers,  walkoff  effects  do  not  occur,  and  therefore  fibers  have  an  additional 
advantage  over  bulk  interactions,  although  fibers  whose  axes  are  not  aligned 
with  the  crystal  axes  may  suffer  leakage  losses. 

Clearly,  single  crystal  fibers  hold  great  promise  for  efficient  low  power 


nonlinear  interactions.  The  discussion  to  this  point  has  been  confined  to  perfect 


waveguides  with  no  axial  variations.  We  will  at  this  time  briefly  discuss  these 
variations  and  the  limitations  they  impose  on  nonlinear  interactions  and  on 
effects  of  ferro  electric  domains. 


*■»  f  u  ^ 

h  uic 


a  Phasematching  Effects 

The  efficiency  of  nonlinear  interactions  depends  critically  on  the  phase 
mismatch,  Ak.  Deviations  in  Ak  from  its  nominal  value  (  Ak=0  for  a 
phasematched  interaction)  can  be  caused  by  axial  variations  in  both  material 
composition  and  fiber  cross-section.  The  effect  of  Ak  variations  are  now 
analyzed  and  related  to  diameter  control  tolerances. 

The  analysis  begins  wirn  rue  equation  for  me  evolution  of  The  amplitude  of 
tne  second  harmonic  mooe 


0-1 '07  =  k  pxpi-iAkZ'  2  i 

wnere  k  is  tne  ooupiino  constant,  and  a(z>  is  the  mooai  amplitude  we  now 
assume  an  axial ly  varying  Ak-  Ak(z)  The  formal  solution  is  straightforward 

31  7  1  =  V.  |pyp( -| /v  |<  1  7  )7' if)?'  ;  | 
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where  we  have  assumed  that  a(0)-0,  i.e.,  that  there  is  no  second  harmonic  input 
at  the  end  of  trio  fiber. 

If  Ak  is  sufficiently  large  that  the  accumulated  phase  difference  becomes 
larger  than  n,  the  second  harmonic  amplitude  is  not  a  monotonically  increasing 
function  of  z,  and  the  efficiency  of  the  interaction  is  severely  degraded.  In  the 
other  limit  where  Ak  is  small,  the  accumulated  phase  difference  remains  small 
compared  to  n,  and  the  second  harmonic  amplitude  increases  monotonically,  Put 
more  slowly  than  for  a  perfectly  phasematched  interaction. 

If  the  accumulated  phase  shift  is  small  compared  to  tt,  we  can  apply  a 
perturbation  approach  and  expand  the  exponential  in  a  power  series,  so  that. 

5(L)  =  K[L  +  i 1 1 (L)  -  i/2l2(L)j  2.13 

where  L  is  the  length  of  the  fiber,  and; 

l  )(L)  =  j  Ak(z')z'dz‘  2M 

b(L)  =  l[Ak(z')z']2dz'  2.15 

Assuming  that  Ak(z)  is  random  with  zero  mean  and  an  autocorrelation  length 
short  compared  to  L,  we  can  integrate  by  parts  equations  2. 1 4  and  2. 1 5  to 
obtain: 


;(L)  =  Lz<Ak>/2  =  0 

2.16 

oil)  =  L^<Ak2>/3 

2.17 

where  Eq.  2  16  follows  from  the  assumption  that  Ak  has  a  zero  mean.  We  can 
now'  see  that  the  second  harmonic  power  P2  is  given  by: 

Poa  ai'D2  =  K2L2[  1  -  (<AK2>/3  )L21  2  18 

We  now  see  that  the  second  harmonic  power  is  reduced  by  (<Ak2>/3)L2  with 
respect  to  a  perfectly  phi  ematched  Interaction.  Therefore,  if  we  set  our 
quality  criterion  by  the  requirement  that  the  power  decrease  be  no  more  than 
one-half  with  respect  to  a  perfect  Interaction,  then  we  must  have. 

<Ak2>  <  (3/2)L"2  2,9 

We  can  relate  this  condition  to  one  on  the  fiber  diameter.  Let  Sp  be  the  RMS 
deviation  of  the  fiber  diameter  from  its  mean  value,  i.e.,  6p=V <p2>,  Then  for 
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small  8p,  we  have; 

5p  <  -HZ/ 2)  (LdAk/dp)- '  2  20 

f- or  fundamental  mode  interactions  in  far  from  cutoff  step  profile  fibers,  we 
have  that 

8p/p  <  (2n/V3)nQ(jQ] )  ^  p^/ (LX)  2.21 

where  jrj i  is  given  by  jo i  =  3tt/4 

For  second  harmonic  generation  with  a  fundamental  wavelength  \  =  1  pm  in  a 
10  crn-lonq,  25  urn  radius  fiber,  we  find  8p/p  <  0.1%  and  it  is  <7%  for  a  1 
em-long  fiber. 

In  the  large  accumulated  phase  difference  limit,  the  decrease  in  efficiency 
depends  in  detail  on  the  form  of  Ak(z),  As  one  might  expect,  the  second 
harmonic  power  is  reduced  by  a  factor  on  the  order  m~2  where  the  accumulated 
phase  difference  is  ~mn. 

b.  Guided  Mode-Guided  Mode  Coupling 

The  modes  in  a  fiber  in  general  have  been  determined  under  the  assumption 
that  the  waveguide  was  2- invar  iant.  In  fibers  with  axial  var  iatioris  in  e,  i.e., 
where  e(x,y,z),  these  modes  no  longer  satisfy  Maxwell's  equations,  ana  corrected 
solutions  must  he  found  The  usual  approach  to  this  problem  is  to  find  solutions 
in  terms  of  a  sum  over  a  complete  set  of  modes,  where  z-dependent  amplitude 
coefficients  represent  the  effects  of  the  axial  variations  in  t  To  simplify  this 
problem,  we  assume  that  the  waveguide  is  isotropic,  so  that  e=n2l. 

The  concept  of  mode  coupling  is  covered  in  Ref,  7.  The  approach  of  local 
mode  expansion  is  used  to  solve  this  problem.  At  each  z.  the  local  mode  is  a 
solution  to  a  hypothetical  z-independent  waveguide  whose  transverse  refractive 
Index  profile  is  the  same  as  that  of  the  perturbed  waveguide  at  that  value  of  z. 
These  z-dependent  modes  are  not  solutions  to  Maxwell's  equations,  and 
therefore  the  true  solutions  must  be  expressed  again  as  a  sum  over  the  local 
modes  with  z-dependent  amplitude  coefficients  that  reflect  the  effects  of  the 
perturbation.  The  advantage  of  the  local  normal  modes  for  problems  involving 
laroe  perturbations  is  that  they  are  discontinuous  at  the  same  radius  as  the 
actual  solution 

The  theory  of  coupled  local  modes  is  developed  in  Ref.  7  and  it  Is  shown  there 
that  Maxwell's  equations  m  the  perturbed  waveguide  can  be  rewritten  in  the 
form  of  a  set  of  coupled  equations  for  the  modal  amplitudes.  The  total  electric 
field  is  expressed  by  the  expansion: 

Et(x.v,z)  =  lCp(z)et(x.y)exp(ik0Z)  2.22 
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The  evolution  of  the  modal  amplitudes  is  described  by. 

dC|j/dz  - 1  Rp/3(z)Q(z)C|3expIi(kj3  -  kjJL)z]  2.23 

We  may  take  k  equal  to  the  conventional  definition  of  k.  represents  the 

strength  of  the  coupling,  while  g(z)  is  related  to  f(z),  the  axial  dependence  of 
the  perturbation  by: 

g(z)  =  -idf/dz  [1/(kfi-kp)]  2.24 

If  we  launch  a  single  mode,  say  13,  into  a  fiber  of  length  L  and  we  further 
assume  that  the  perturbations  are  sufficiently  weak  that  only  a  small  fraction 
of  the  incident  power  is  scattered  into  other  modes,  then  we  can  assume  that  to 
first  order,  C^L)  =  C^tO).  Eq.  2.23  can  be  easily  integrated  to  find  the 

coefficients  of  the  other  modes: 

Cjjt'L)  =  R^p,  Cp(0)  \  g(z)  exp[i(kj3  -  k^zldz  2.25 

integrating  this  equation  once  by  parts  and  assuming  that  the  perturbation 
vanishes  at  0  and  L,  we  find: 

Cjjd)  =  RuB  Cb(0)A  F(kB-k(J)  2.26 

where  F(k)  is  the  Fourier  transform  of  the  perturbation  function  f(z) 

F(k)  =  ( l/VL)  J  f(z)  exp(ikz)dz  2.27 

We  then  find  that  the  ratio  of  the  power  carried  by  each  mode  is  given  by. 
Pm(L)/Pb(0)  =  IRui3t2F(kB-ku)2L  2.28 


Clearlv,  strong  coupling  occurs  only  when  the  Fourier  spectrum  of  the 
perturbation  contains  components  at  the  difference  frequency  of  the  coupled 
modes.  We  can  express  a  periodic  perturbation  as  having  the  form: 

f(z)  =  sin(2nz/A)  2.29 

We  define  a  quality  criteria  for  perturbations  as  the  loss  over  a  length  L  of 
less  tnan  a  fraction  t  or  the  power  incident  on  the  waveguioe  For  a  resonant 
p^nodK  perturbation,  we  have  that: 
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IR^I2  <  4t/L2 


2.30 


if  we  assume  the  coupling  to  be  due  to  pure  diameter  variations  on  a  step 
profile  fiber  and  that  the  mode  is  far  from  cutoff,  we  can  calculate  the 
tolerance  for  coupling  from  an  HE]  ]  mode  to  an  HE  j , n  mode  3S: 

6p/p  <  f  1.7/(m-0.25)][nc/'h(A)](  p2/LX)Ve  2.31 

For  X=  1  pm,  p=  1 0pm,  nc=2  and  L=  1  cm,  we  have  that: 


ftp/ p  <  0.034  Ve/[(m-0.25)h(A)]  2.32 

where  h(A)  varies  between  1  for  weakly  guiding  fibers  to  roughly  5  for  unclad 
fibers.  If  we  take  e=0.0 1 , 1.e.,  1%  loss  in  I  cm,  we  find  that  for  clad  fibers,  the 
diameter  variations  must  be  smaller  than  several  tenths  of  a  percent  for 
coupling  tc  the  HE] 2  mode,  and  smaller  roughly  as  1/m  for  coupling  to  higher 

order  HE  ] m  modes  it  is  also  useful  to  know  the  scale  of  perturbation  that 
causes  coupling  from  the  HE  ]  |  mode  to  the  HE ] m  mode.  The  resonance  condition 
for  far  from  the  cutoff  modes  is  given  by 

a  =  8nc/l(m-0  25 >2  -  (0.75 )2](  p2/X)  2  33 

With  the  same  numerical  parameters  used  in  Eq.  2.32,  we  find. 

A(prn)  =  1 600/[(rn-0.25)2  -  (0.75)2]  2.34 

For  rne  he  ]2  mode,  we  find  A  =  600pm,  and  tne  resonant  wavelength  decreases 
roughly  as  1  /m2  for  coupling  to  higher  order  HE  ]  m  modes. 

wmie  this  modei  is  crude,  it  can  serve  as  an  aid  in  correlating  experimental 
data  We  note  in  passing  that  X"1  loss  scaling  was,  in  fact,  observed  in 
measurements  of  radiation  losses  of  sapphire  fibers.  At  Stanford,  It  was  also 
observed  that  losses  in  diameter  controlled  fibers  were  not  improved  over  those 
in  good  quality  open  loop  fibers  of  the  same  diameter,  despite  a  several-fold 
reduction  in  the  RMS  deviation.  This  unexpected  situation  can  be  explained  by 
the  increased  spatial  frequency  observed  for  the  perturbations  in  the  closed 
loop  fibers,  and  the  A-2  scaling  of  the  loss.  To  completely  explain  this 
question,  quantitative  power  spectra  must  be  determined,  but  the  implications 
of  this  simple  model  can  be  useful  in  interpreting  any  data. 


c.  Radiation  Losses 


The  other  important  loss  mechanism  caused  by  structural  perturbations  is 
the  coupling  from  the  desired  guided  mode  to  the  radiation  modes.  Calculation 
of  these  radiation  losses  is  much  more  difficult  than  for  guided  mode  losses. 
The  most  common  approach  is  a  coupled  mode  analysis.  For  weakly  guiding 
fibers,  the  coupling  coefficients  are  considerably  more  complicated  than  those 
for  guided  modes.  For  unclad  fibers,  one  must  calculate  coupling  coefficients 
from  the  exact  radiation  modes  of  a  cylinder  with  arbitrary  discontinuities  in 
the  dielectric  constant.  With  the  exception  of  the  TE  modes,  this  analysis  is 
extremely  complex  and  has  been  accomplished  only  for  these  and  the  HE  j  \ 

modes8-10.  For  large  core-cladding  index  differences,  other  approximations 
have  been  applied,  including  ray  tracing  approaches1 J'12,  Fokker-Planck 
diffusion-like  equations 4  and  analogies  to  bulk  scattering1 5'  17  induced 
current  representations12  have  also  been  used.  Finally,  radiation  losses  caused 
by  perturbations  with  too  narrow  a  power  spectrum  to  couple  directly  to 
radiation  modes  have  been  analyzed  in  terms  of  the  statistical  probability  that 
a  particular  guided  mode  fal  Is  into  the  cutoff  regime. 1 8 

We  confine  our  discussion  to  qualitative  comments  based  on  the  literature 
results  for  radiation  losses.  From  the  synchronicity  condition,  it  is  clear  that 
radiation  mod*3  coupling  is  mediated  by  components  of  the  perturbation  that  are 
shorter  than  those  causing  guided  mode  coupling.  It  has  been  shown  that  in 
planar  waveguides  with  random,  exponentially  autocorrelated  perturbations,  for 
both  large  and  small  Index  differences,  when  fewer  than  three  modes  are  below 
cutoff,  the  peak,  loss  of  the  fundamental  mode  is  dominated  by  radiation  losses, 
while  for  three  guided  modes  the  peak,  loss  is  dominated  by  coupling  to  guided 
modes.  Since  the  loss  results  are  generally  similar  for  fibers  and  slab 
waveguides7'0,  it  is  reasonable  to  expect  that  losses  for  multimoded  f  ibers  are 
dominated  by  modal  scattering  rather  than  radiation  losses.  Thus,  we  may 
assume  that  the  diameter  tolerances  derived  assuming  only  guided  mode  losses 
are  adequate,  unless  the  perturbation  power  spectrum  is  very  wide,  which  has 
not  been  observed 

d.  Ferroelectric  Domains 

e. s  is  typical  of  ferroelectric  materials,  as-grown  crystals  tend  to  break  up 
into  a  number  of  regions  of  reversed  spontaneous  polarization  < ferroelectric 
domains)  Since  the  sign  of  most  of  the  material  parameters  differs  in 
antiparailel  domains,  the  efficiency  of  devices  depends  on  the  preparation  of 
"poled"  fibers 

As-grown  bulk  crystals  of  LiNbCH  tend  to  exhibit  multi-domain  structures 

These  domain  distributions  are  conventionally  explained  in  terms  of  a  balance 
depolarization  energy,  which  favors  a  single  domain  crystal4  In  c-axis  rods, 
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between  these  domains  tend  to  take  the  form  of  concentric  rings,  while  in 
a-axis  rods  the  domains  form  a  series  of  layers  perpendicular  to  the  rod  axis. 
The  domain  wall  spacing  that  minimizes  the  sum  of  the  depolarization  energy 
and  the  domain  wall  energy  in  an  a-axis  rod  has  been  calculated  by  M.  Fejer.  The 
geometry  considered  is  shown  in  Fig.  2.1.  The  result  for  the  domain  wall 
spacing.  A,  as  a  function  of  the  rod  radius,  a,  is  found  to  be. 

A  =  Vra((  TT3/2  10)e.o(f  ♦  1  )c/Ps2]1/2  2  3b 

where  t  is  the  dielectric  constant  of  the  rod,  c  is  the  domain  wall  energy,  and  Ps 

is  the  sponteneous  polarization  per  unit  volume.  Since  these  material 
properties  are  not  well  known  near  the  Curie  temperature,  we  can  best  estimate 
the  predicted  domain  wall  spacing  in  fibers  by  noting  that  the  quantity  in  square 
brackets  is  independent  of  rod  diameter,  and  simply  scaling  the  known  spacing 
in  bulk  crystals  by  the  square  root  of  the  diameter.  Typical  domain  size  for 
centimeter  diameter  bulk  a-axis  boules  is  several  hundred  microns.  We  might 
thus  expect  domain  spacing  of  50-  100pm  in  500pm  diameter  fibers. 

Fejer  measured  the  domain  distribution  in  L1Nb03  fibers  by  measuring  the 

change  in  the  sign  of  the  pyroelectric  coefficient  in  antiparallel  domains.  The 
voltage  developed  between  a  fine  wire  on  one  side  of  the  crystal  and  a  ground 
plane  on  the  other  side  of  the  crystal  was  monitored  while  the  surface  of  the 
fiber  was  heated  with  a  chopped  CO2  laser  beam.  The  relative  phase  between 

the  chopper  and  the  voltage  across  the  crystal  was  taken  as  an  indication  of  the 
sign  of  the  pyroelectric  coefficient,  and  thus  of  the  orientation  of  the  domain. 
The  distribution  of  domains  along  the  axis  could  then  be  measured  by  sliding  the 
wire  along  the  fiber  This  technique  clearly  Is  most  applicable  to  fibers  whose 
axis  is  normal  to  the  c-asix,  and  depends  on  an  assumption  that  the  domains  are 
uniform  along  the  c-axis.  This  assumption  seemed  reasonable  in  terms  of  the 
energy  argument  presented  above,  as  head-to-head  domains  are  a  mgher  energy 
configuration  than  no  domains  at  all  They  were  surprised  when  the  phase  of  the 
voltage  produced  by  an  a-axis  fiber  was  found  to  be  independent  of  the  position 
along  the  axis,  and  invariant  under  a  180'  rotation  of  the  fiber.  The  latter 
observation  suggested  that  head-to-head  domains  were  forming  along  the 
c-axis,  while  the  former  suggested  that  the  fiber  did  not  break  up  into 
thelayered  domain  structure  characteristic  of  bulk,  crystals.  The  results  for  the 
fibers  to  be  essentially  similar  to  those  in  the  bulk,  they  at  first  attributed 
these  results  to  experimental  error,  However,  these  measurements  were 
verified  both  with  a  piezoelectric  technique,  and  a  more  sophisticated 
pyroelectric  apparatus19,  The  pyroelectric  device  Is  capable  of  mapping  out  the 
domain  distribution  across  the  face  of  a  fiber.  They  found  that  a-axis  fibers 
did,  in  fact,  break  up  into  head-to-head  domains,  with  the  domain  wall  roughly 
at  the  fiber  axis  These  results  for  a-axis  fibers  were  further  verified  by 


etching  studies  by  Nightingale.20  Luh  made  a  series  of  careful  etching  studies 
of  both  a-  and  c-axis  fibers.21  He  again  found  the  same  results  for  a-axls 
fibers,  and  discovered  that  c-axis  fibers  up  to  800|inn  In  diameter  were  single 
domain.  The  positive  c-axis  of  the  fiber  was  always  found  to  point  towards  the 
melt,  independent  of  seed  orientation. 

It  is  clear  that  simple  energy  considerations,  including  only  depolarization 
and  domain  wall  energy,  are  unable  to  account  for  these  empirical  results. 

Other  energy  terms  that  affect  the  domain  distribution  must  be  considered  it 
is  known  that  the  spontaneous  polarization  in  ferroelectrics  tends  to  align  with 
externally  applied  electric  fields.  At  temperatures  near  the  melting  point, 
fields  as  small  as  0.4V/cm  applied  along  the  c-axis  of  bulk  crystals  have  been 
found  to  yield  a  poled,  single  domain  structure.22  Thus,  the  domains  observed  in 
crystal  fibers  could  be  explained  by  an  internally  generated  electric  field.  A 
mechanism  for  the  generation  of  such  fields,  the  thermoelectric  effect,  exists 
and  can  serve  as  the  basis  for  a  model  that  explains  all  observations  as  shown 
by  Fejer. 

The  electric  field  generated  by  carriers  diffusing  in  the  influence  of  a 
thermal  gradient  can  be  written  in  the  form; 

E  =  Q  ■  VT  2.36 

where  E  is  the  electric  field,  T  Is  the  tempeature,  and  Q  Is  the  second  rank 
thermoelectric  power  tensor.  For  our  purposes,  it  Is  adequate  to  assume  that  Q 
is  isotropic  and  of  magnitude  Q.  Measurement  of  Q  in  a  bulk  sample  of  congruent 
LiNbOj  in  air  at  temperatures  close  to  the  Curie  temperature  yielded  a  value  of 

O.emV/K21 

Consider  now  the  growth  of  a  c-axis  fiber.  For  our  condition  of  fiber  growth, 
the  temperature  gradient  along  the  axis  of  the  fiber  can  be  estimated 
conservatively  as  greater  than  1000  K/cm.  Using  Eq.  2.36,  we  can  determine  an 
axial  electric  field  larger  than  1  V/cm,  greater  than  that  needed  to  pole  bulk 
crystals.  Further,  the  orientation  of  the  electric  field  in  independent  of  the 
orientationof  the  seed  crystal,  which  explains  the  preferred  orientation  of  the 
domains  in  c-axis  fibers  In  a-axis  fibers,  the  radial  temperature  gradients  lead 
to  a  radial  electric  field.  The  projection  of  the  electric  field  on  the  c-axis  is 
then  of  precisely  the  form  necessary  to  explain  the  head-to-head  domains 
observed  in  a-axls  fibers. 

The  model  was  further  tested  by  heating  an  a-axis  fiber  asymmetrically  to 
shift  the  temperature  gradient  from  a  symmetric  to  an  asymmetric  position 
v.'ith  respect  to  the  fiber  axis.  The  position  of  the  domain  wall  w*s  feund  to 
move  away  from  the  center  of  the  fiber,  as  predicted  by  the  thermoelectric 
model  21 

As  pointed  out  previously,  it  is  desirable  from  a  device  standpoint  to  use 
single-domain  fibers,  As  mentioned,  the  as-grown  c-axis  fibers  fortuitously 


take  this  form.  The  results  of  the  asymmetric  heating  experiments  just 
described  suggest  an  approach  to  poling  an  a-axls  fiber.  It  Is  clear  from  the 
thermoelectric  poling  model  that  the  direction  of  the  spontaneous  polarization 
in  a-axis  fibers  is  dictated  by  the  sign  of  the  transverse  temperature  gradient. 
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Section  3 


LASER  HEATED  PEDESTAL  GROWTH  TECHNIQUE 

It  is  clear  from  the  discussion  given  earlier  that  the  growth  of  crystal 
fibers  of  high  melting  point  materials  must  overcome  several  difficult 
challenges  in  order  to  be  useful  for  the  many  applications  listed.  A  number  of 
techniques  have  been  used  to  produce  single  crystal  fibers  and  these  are 
schematically  illustrated  in  Figure  3.1. 

In  the  hot  rolling  process,  shown  in  Figure  3. 1  (a),  a  rod  is  extruded 
through  a  series  of  progressively  smaller  wire  dies  to  produce  a  fiber  of 
reduced  diameter.25  The  method  is  primarily  applicable  to  relatively  soft 
materials  with  low  melting  points,  e.g.  alkali  halides. 

In  the  capi  1  lary-Br idgeman  approach,  Figure  3.1(b),  a  glass  capillary  tube  is 
immersed  in  a  melt  of  the  crystal  to  be  grown.  The  apparatus  is  slowly 
lowered  through  a  temperature  gradient  to  produce  a  crystalline  core  inside 
the  glass  tube24  Capillary-Bridgeman  growth  has  primarily  been  applied  to 
low-melting  point  organic  crystals,  as  thermal  expansion  matching  and 
chemical  compatibility  are  critical  requirements  in  this  technique. 

A  number  of  techniques,  for  example  the  edge-defined  film-fed  growth  and 
related  processes,  involve  feeding  the  melt  through  a  capillary  to  a  die,  which 
serves  to  shape  the  resulting  fiber25,  Figure  3.1(c).  These  methods  are 
reviewed  in  reference  26.  Single  crystal  fibers  of  silver  and  thallium  halides 
have  been  produced  by  such  techniques27'28,  as  have  lower 
quality  fibers  of  high  melting  point  oxides  such  as  sapphire,  lithium  niobate 
and  spinel  25  Particularly  for  high  melting  point  materials,  the  choice  of 
appropriate  die  material  is  crucial  to  minimize  contamination  of  the  melt  and 
provide  proper  wetting  conditions  for  stable  growth.  For  sapphire 
growth  with  such  a  high  melting  point,  this  technique  leads  to  in-diffusion  of 
the  die  metal  such  as  molybdenum. 

Laser-heated  pedestal  growth  (LHPG),  illustrated  in  FI  jure  3.1(d),  uses  a 
focussed  laser  beam  to  melt  the  tip  of  a  solid  rod  of  the  material  to  be  grown. 
A  seed  crystal  is  dipped  into  the  melt,  then  withdrawn  at  a  faster  rate  than 
the  source  rod  is  fed  in,  resulting  in  the  growth  of  a  fiber  of  reduced 
diameter.  This  approach  has  the  advantages  of  eliminating  any  possibility  of 
crucible  contamination  and  of  achieving  melt  temperatures  limited  only  by 
the  available  laser  power.  Moreover,  difficulties  in  the  stabilization  of  the 
fiber  diameter  have  in  large  part  been  solved.  LHPG  has  been  applied 
primarily  to  refractory  oxide  materials29-5 1 ,  though  halides,  borides,  and 
elemental  semiconductors  have  also  been  grown  52 

No  one  of  these  approaches  is  best  suited  to  the  production  of  high  quality 
fibers  in  all  materials.  We  have  focused  on  the  refractory  oxide  materials, 
for  which  LHPG  offers  great  flexibility  and  has  the  potential  to  produce  the 
purest  fibers. 


17 


Figure  3.1.  Techniques  for  the  growth  of  single  crystal  fibers,  (a)  hot 
rolling;  b)  Bridgeman  growth  in  glass  capillaries; 

(c)  edge-defined  film-fed  growth;  (d)  laser-heated 
pedestal-growth  (LHPG) 
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As  mentioned  earlier,  the  LHPG  technique  involves  melting  the  tip  of  a 
source  rod  with  a  focussed  laser  beam,  dipping  in  an  oriented  seed  crystal, 
then  pulling  a  fiber  while  feeding  in  the  source  rod.  This  process,  quite 
simple  in  principle,  has  several  attractive  features  from  a  materials 
processing  standpoint. 

( 1 )  Very  high  temperatures  can  be  achieved,  limited  in  principle  only  by 
the  radiation  temperature  of  the  laser,  and  in  practice  by  the  available  laser 
power  Only  a  few  watts  of  absorbed  laser  power  are  necessary  to  raise  the 
temperature  of  500  pm  diameter  rods  of  typical  oxide  materials  to  their 
melting  point. 

(2)  The  temperature  gradients  at  the  freezing  interface  are  on  the  order  of 
thousands  of  degrees  per  cm.  These  large  temperature  gradients,  in  turn, 
allow  stable  growth  at  rates  orders  of  magnitude  larger  than  in  bulk  crystals, 
typically  0. 1  -20mm/min.  As  thermal  stresses  in  a  surface  cooled  cylinder 
scale  with  the  radius,  these  large  gradients  do  not  lead  to  a  cracking  problem 
in  small  diameter  fibers. 

(3)  The  process  is  entirely  crucibleless,  minimizing  the  incorporation  of 
undesired  dopants,  and  facilitating  the  growth  of  a  variety  of  materials  in  the 
same  apparatus  with  a  high  degree  of  purity  equal  to  the  purity  of  the 
starting  material. 

(4)  It  is  theoretically  predicted  that  in  the  absence  of  phase  separation  or 
volatilization,  the  composition  of  the  source  rod  and  the  fiber  should  be 
identical,  so  that  the  effective  distribution  coefficient  approaches 

unity-  —  4,  which  is  particularly  important  in  the  growth  of  compositional ly 
uniform  solid-solution  crystals.  The  prediction  that  convection  in  the  molten 
zone  is  not  oscillatory  further  reduces  the  probability  of  compositional 
striations 

In  order  to  take  advantage  of  these  features  of  the  LHPG  process  to 
produce  fibers  of  high  melting  point  materials,  stable  growth  conditions  must 
be  realized,  which  in  turn  requires  thermal  and  mechanical  stability  of  the 
growth  zone,  and  symmetric  heat  input  from  the  laser.  The  heat  should  be 
delivered  to  a  spot  whose  size  is  comparable  to  the  source  rod  diameter,  to 
produce  the  short  molten  zones  necessary  for  stable  growth.  The  apparatus 
designed  to  meet  these  requirements  is  now  described. 

A  block  diagram  of  our  fiber  growth  apparatus  is  shown  in  Figure  3.2.  A 
focused  CCb  laser  melts  a  surface  tension  supported  liquid  zone  which 

bridges  the  source  ana  seed  rods  Growth  proceeds  by  simultaneous  upward 
translation  of  the  seed  and  source  rods  with  the  molten  zone  positioned 
between  them  The  laser  focal  spot,  and  consequently  the  molten  zone,  remain 
fixed  during  fiber  growth  The  source  rod  to  fiber  diameter  ratio  is  set  by 
mass  conservation  to  be  the  square  root  of  the  fiber  to  source  rod  translation 
rate 

in  order  to  achieve  a  constant  fiber  diameter,  stable  fiber  growth 
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conditions  must  be  realized  This  in  turn  dictates  a  rigid  mechanical 
apparatus,  smooth  source  feed  and  fiber  pull  rates,  stable  laser  power,  and 
symmetric  heat  input  into  the  molten  zone.  The  LHPG  apparatus  utilizes  novel 
optical  and  mechanical  systems  to  achieve  stable  growth  conditions  and  a 
uniform  diameter  fiber.  Descriptions  of  the  fiber  growth  apparatus 
sub-systems  are  given  below. 

A  polarized  CO9  laser  serves  as  the  heat  source  for  crystal  growth.  The 

water  cooled  laser  cavity  is  temperature  stabilized  and  produces  a  polarized 
HF  ]  |  output  mode  with  power  fluctuations  of  less  than  0.75%.  A  polarization 

power  control  system  is  used  to  adjust  the  laser  power  incident  on  the  molten 
zone.  After  passing  through  a  ZnSe  beam  expanding  telescope  and 
some  beam  steering  optics  the  CO2  beam  enters  the  controlled  atmosphere 
orowth  chamber. 

within  the  growth  chamber  a  novel  optical  system  focuses  the  laser  beam 
onto  tr,e  fiber  with  3  360  degree  axially  symmetric  distribution  as  shown  in 
Figure  3  3  The  symmetric  irradiance  prevents  cold  spots  in  the  growth  zone 
ano  represents  a  significant  improvement  over  the  previously  used  two 
beam  35  rotating  periscope, 36  or  ellipsoidal-"  focusing  systems 

a  novel  optical  element  incorporated  into  the  design  is  a  reflexicon^ 
which  rpnsists  of  an  inner  cone  surrounded  by  a  larger  coaxial  cone  in  order 
to  acmeve  good  optical  performance  it  is  critical  that  the  reflexicon's  two 
cones  be  accurately  aligned  A  mated  surface  design  using  diamond  turned 
copper  optical  components,  assures  centering  of  the  cone's  axes  A  gold 
coating  on  the  copper  optical  surfaces  enhances  reflectivity  ana  protects  the 
copper  substrate  The  ref  lexicon  and  parabolic  mirror  provide  near  diffraction 
limited  f/2  focusing,  yielding  a  minimum  spot  size  of  30  microns  This  tight 
focus  is  important  for  the  stable  growth  of  small  diameter  fibers  The  focal 
spot  size  can  be  controlled  by  modifying  the  input  beam  divergence  with  the 
focusing  telescope  X-Y  stages  on  the  fiber  and  source  rod  translation  devices 
permit  adjustment  of  the  fiber  position  with  respect  to  the  fixed  laser  focal 
spot 

Fiber  translation  speed  is  controlled  by  a  phase-locked  control  circuit 
which  enables  d  c  motor  operation  over  a  250,000  1  speed  range  with 
acquisition  times  of  1  ms  A  useful  control  option  allows  the  pull-to-feed 
translation  ratio  to  be  fixed  while  adjusting  the  growth  speed  Since  starting 
transients  evjst  during  initiation  of  fiber  growth,  a  slow  growth  ratp  allows 
time  for  adjustments  As  the  equilibrium  growth  conditions  are  reached,  the 
growth  <;pppg  can  be  increased  without  affecting  the  fiber  diameter 

The  first  step  in  the  growth  of  a  fiber  is  the  preparation  of  the  source 
material  The  requirements  on  the  source  rod  are  (Pit  should  have  the 
composition  of  the  oesirer)  end  product  (assuming  that  no  volatilization  fakes 
place  during  growth',  (?)  it  should  have  constant  density,  and  have  as  close  to 
theoretical  density  as  possible,  (3)  the  cross-section  should  be  constant,  as 
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Figure  3.2.  Block  diagram  of  the  LHPG  apparatus. 
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Figure  3.3.  Cross-section  of  the  growth  head,  showing  the  reflexicon 
focussing  system  and  the  fiber  translation  mechanisms.  The 
optical  elements  are  rotationally  symmetric  about  their 
centerline. 


small  as  possible,  and  preferably  circular.  We  have  had  good  success  using  a 
centerless  grinder  to  fabricate  starting  rods  with  diameters  as  small  as 
300  pm  and  tapers  as  small  as  1  pm/cm.  The  rod  material  itself  can  be  a 
single  crystal  obtained  elsewhere,  a  polycrystalline  source  material  from 
solidified  melts,  a  cold-pressed  and  sintered  powder,  or  a  hot-pressed  powder. 
All  of  these  have  been  used  successfully  20,39,40 

The  next  step  in  the  growth  is  to  melt  the  tip  of  the  source  rod  with  the 
CO>  laser  beam,  and  to  dip  in  a  seed  crystal.  The  seed  can  be  an  oriented 

centerless-ground  rod,  a  previously  grown  fiber  or,  in  the  case  of  a  first 
growth,  a  platinum  wire  or  a  fiber  of  a  higher  melting  point  material.  Because 
of  the  symmetrical  heat  input  of  the  reflexicon  focussing  system,  the 
azimuthal  orientation  of  the  seed  has  no  direct  effect  on  the  growth. 

After  the  seed  is  oriented,  the  laser  power  is  adjusted  to  produce  the 
desired  length  of  molten  zone.  The  molten  zone  must  also  be  adjusted  to 
produce  the  proper  contact  angle  between  the  melt  and  the  seed 

Growth  is  initiated  at  this  point  by  simply  simultaneously  switching  on  the 
pull  and  feed  motors,  which  have  previously  been  set  to  run  at  the  appropriate 
speeds  Lengths  over  100  cm  have  been  grown  to  date.  In  principle,  lengths  of 
fiber  are  limited  only  by  the  available  feed  material.  This  latter  limitation  is 
not  serious  in  practice,  in  that  each  centimeter  length  of  1 .25  mm  diameter 
feed  material  could  produce  5  meters  of  50  jim  diameter  fiber  if  processed 
three  times 

Only  a  modest  amount  of  power  is  necessary  to  produce  an  appropriate 
molten  zone  For  example,  2  watts  incident  power  is  adequate  for  the  growth 
of  a  1 70|im  diameter  sapphire  fiber  from  a  500pm  diameter  source  rod 
typical  growth  rates  are  0.1-20  mm/min,  several  orders  of  magnitude  faster 
than  is  typical  of  bulk  crystal  growth  The  mechanism  that  limits  the 
maximum  possible  growth  rate  appears  to  vary  for  different  materials  in 
Reference  20,  Nightingale  shows  that  for  sapphire  grown  at  rates  faster  than 
8  mm/min,  constitutional  supercooling  with  respect  to  an  unknown  species, 
possibly  At  or  O2  leads  to  the  formation  of  microvoids  along  the  axis  of  the 

fiber  that  cause  severe  optical  scattering  problems.  Growth  in  a  He 
atmosphere  to  increase  the  thermal  gradients  at  the  Interface  may  lead  to 
higher  useful  growth  rates. 

We  have  generally  found  it  possible  to  grow  fibers  approximately  a  factor 
of  three  smaller  in  diameter  than  the  source  rod  for  growth  in  air.  As  the 
diameter  reduction  is  increased,  the  damping  coefficient  for  diameter 
variation  decreases,  eventually  becoming  negative  at  a  critical  diameter 
reduction  Rr.  Growth  at  diameter  reductions  larger  than  Rc  is  thus  unstable. 

For  AI2O3  this  instability  occurs  at  a  reduction  between  3  and  4  when  grown 

in  air  it  is  possible  that  diameter  reductions  between  A  and  5  could  be 
achieved  when  qrown  in  a  helium  atmosphere 

LaserGenics  Corporation  has  recently  obtained  the  exclusive  rights  to  the 
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three  patents  held  by  Stanford  University  on  this  single  crystal  fiber  growth 
process  using  the  LHPG  technique.  We  are  developing  this  process  in  order  to 
obtain  high  quality  single  crystal  fibers  of  several  meter  lengths.  The  LHPG 
technique  is  an  excellent  method  to  grow  single  crystal  fibers.  It  Is  likely 
that  for  many  high  melting  point  materials  no  other  technique  will  be  able  to 
grow  similarly  good  quality  single  crystal  fibers. 

With  this  process,  the  atmosphere  can  be  accurately  controlled  during 
growth  as  well  as  the  crystal  growth  temperature.  Because  of  the  steep 
thermal  gradients  In  the  growth  zone  with  this  technique,  the  crystal 
structure  that  Is  "frozen  In"  can  to  some  degree  be  controlled.  We  feel  that 
this  Is  an  Ideal  technique  to  Investigate  the  growth  of  high  melting  point 
materials.  This  technique  is  closely  related  to  the  float-zone  method  which  Is 
the  best  known  technique  to  grow  incongruently  melting  compositions. 
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SECTION  A 


RESULTS  AND  CONCLUSIONS 

In  this  section  we  will  first  review  the  results  that  have  been  obtained  at 
Stanford  University.  The  purpose  is  to  inform  the  reader  of  the  earlier  results 
and  to  demonstrate  that  the  theory  described  earlier  is  indeed  valid.  These 
results  have  been  obtained  primarily  by  Dr.  Marty  Fejer,  Dr.  John  Nightingale  and 
others  of  the  Stanford  group.  Most  of  these  results  have  been  obtained  on 
sapphire  and  YAG  fibers,  however,  they  do  validate  the  theory  presented. 
Following  this  discription  we  will  present  our  results  on  the  growth  of 
MgO.LiNbOj  fibers. 

As  remarked  earlier,  the  optical  performance  of  devices  made  of 
single-crystal  fibers  is  controlled  in  large  part  by  the  scatter  losses.  As  we 
noted  earlier,  three  types  of  losses  are  important  In  single-crystal  fibers: 
guided  mode  coupling,  radiation  losses,  and  material  absorption.  The  latter  two 
are  readily  measured  by  integrating  sphere  and  calorimetric  techniques, 
respectively,  but  the  guided  mode  coupling  is  far  more  difficult  to  quantify 
Unfortunately,  the  guided  mode  coupling  is  also  expected  to  be  the  largest 
component  of  the  total  losses.  It  is  therefore  important  to  obtain  at  least  an 
estimate  of  the  guided  mode  loss. 

a  Guided-Mode  Coupling 

Tne  sketchy  available  data  are  the  result  of  measurements  made  by  Digonnet 
and  Gaeta  on  single-crystal  Nd  YAG  fiber  lasers,  given  in  Ref.  41.  They  made 
lasers  in  both  open-loop  and  closed-loop  fibers,  where  the  closed-loop  fibers 
wem  gmwn  under  feedback  control  and  typically  showed  0.5%  RMS  diameter 
variations,  while  the  open  loop  fibers  were  grown  without  the  benefit  of 
feedback  control,  and  typically  had  2%  RMS  diameter  variations  Quantitative 
power  spectra  of  the  diameter  variations  are  not  available  for  these  fibers,  so 
bnth  these  amplitudes  and  their  estimated  autocorrelation  length  of  10  x  (fiber 
radius)  should  be  regarded  as  approximate.  Both  clad  and  unclad  fibers  were 
used,  where  the  clad  fibers  were  extrusion  coated  with  LaSF30  glass,  leading  to 
a  cnre-ciadaing  index  of  refraction  difference  of  0.048. 

The  data  can  be  summarized  as  follows:  For  a  series  of  unclad  open  loop 
fibers  with  diameters  ranging  between  75  and  225pm,  the  loss  scaled 
approximately  as  p"3/2(  where  p  is  the  fiber  radius.  If  the  loss  is  expressed  in 
cm'1  and  the  fiber  radius  in  cm,  the  following  expression  is  a  good  fit  to  the 
data: 


a-Cp-3/2 


of 


4.1 


where  C  *  1 .25  x  1 O'4  cm 1  though  there  is  considerable  scatter  in  the  data. 
The  loss  of  an  unclad  closed  loop  fiber  of  48pm  diameter  was  measured  as 
0.06cm"',  while  a  41  urn  diameter  clad  open-loop  fiber  had  0.02cm'1  loss. 

To  compare  these  results  with  theoretical  predictions,  consider  Fig  4. 1 
which  plots  the  normalized  loss  for  unclad  Nd:YAG  fibers  vs.  the  normalized 
autocorrelation  length  D/p,  using  a  profile  height  of  a=0.35  and  an 
autocorrelation  decay  exponent  s=3.  The  precise  choice  of  s  does  not  greatly 
affect  the  results  for  the  cases  discussed  here,  a  Is  defined  as  [  1  = 

(nc]/n0)2]/2  The  loss  a  is  normalized  according  to  ap3/[oh(A)]2,  where  a  is  the 

RMS  diameter  variation  and  h(A)  is  the  correction  factor  for  large  index  of 
refraction  differences,  The  three  curves  represent  different  values  of  2nnp/X, 
where  n  is  the  index  of  refraction  of  the  core,  and  \  Is  the  free  space 
wavelength  For  Nd:YA6  with  n=  1 .8  and  X=  1 .06pm,  the  curves  correspond  to  p 
=25,  50  and  100pm  Two  general  trends  can  be  extracted  from  this  graph  for 
0.02<D/p<20.  First,  the  losses  scale  as  p"3/2,  in  agreement  with  experimental 
observations  for  the  unclad  open-loop  fibers.  Second,  the  loss  depends  only  on 
(D/p) ' /2,  so  that  small  errors  in  estimating  this  parameter  do  not  seriously 
affect  the  results.  Losses  calculated  according  to  the  fit  to  the  data,  Eq.  4. 1 , 
are  plotted  on  Fig  4.1,  and  can  be  seen  to  fall  close  to  the  theoretical  curves. 

The  loss  for  the  unclad  closed-loop  fiber  is  also  plotted  in  Fig.  4. 1 ,  and  is  again 
seen  to  fall  on  the  theoretical  curves. 

Now  consider  Fig  4.2,  which  plots  normalized  loss  calculated  as  a  function 
of  normalized  autocorrelation  length  D/p  for  a  fixed  normalized  radius  of 
2rmp'A  =  267  This  normalized  radius  corresponds  to  50pm  diameter  at  1.06pm 
wavelength  The  two  curves  are  for  different  values  of  the  profile  height  a, 
appropriate  for  unclad  and  cladded  Nd:YAG  fibers.  The  normalized  curves  are 
indistinguishable  over  the  range  of  D/p  of  interest,  so  the  dependence  of  loss  on 
cladding  index  is  only  through  h(A).  The  loss  of  the  clad  41  pm  diameter  fiber 
can  be  compared  to  that  of  the  48pm  unclad  fiber  using  the  c^hCA^p"^2 
scaling  given  previously.  The  loss  of  the  unclad  fiber  is  predicted  to  be  1.5 
times  larger  than  that  of  the  clad  fiber,  compared  to  the  experimentally 
observed  factor  of  3. 

The  above  results  indicate  that  the  simple  correlations  predicted  by  theory 
are  able  to  scale  the  losses  for  a  factor  of  5  variation  in  radius,  a  factor  of  12 
in  index  profile  depth,  and  a  factor  of  5  In  RMS  diameter  variation.  The 
experimental  data  are  uncertain  by  factors  on  the  order  of  two,  so  that  the 
accuracy  of  the  theoretical  predictions  can  be  taken  to  agree  with  experiments 
only  within  this  limit.  Fiowever,  given  the  large  range  of  parameters  over  which 
the  fits  are  reasonable,  the  theory  can  be  regarded  as  correctly  reproducing  the 
essential  dependencies. 
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Figure  4.1.  Mode-coupling  loss  normalized  to  <AAjp3  vs.  the  normalized 
autocorrelation  length.  (Curves  labled  with  the  value  of  2fmp/X) 
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Figure  4.2.  Mode -coupling  loss  normalized  to  (A2/p3  vs.  the  normalized 
autocorrelation  length  at  a  fixed  ratio  of  fiber  radius  to 
wavelength. 


b.  Radiation  Losses 


Guided-mode  coupling,  discussed  in  the  previous  subsection,  is  the  loss 
mechanism  of  primary  importance  in  nonlinear  devices,  where  light  removed 
from  the  fundamental  mode  is  essentially  lost  to  the  interaction.  For  power 
transmission  applications,  one  is  less  interested  in  the  distribution  of  the 
energy  among  the  modes  than  the  power  that  actually  leaks  from  the  fiber. 

Thus,  even  though  the  coupling  to  radiation  modes  is  expected  to  be  smaller  than 
that  to  guided  modes  for  fibers  whose  diameter  variations  have  narrow  power 
spectra,  the  radiation  losses  are  important  in  these  passive  applications. 

The  coupling  between  the  fundamental  mode  and  radiation  modes  is  dictated 
by  the  high  spatial  freguency  components  of  the  diameter  variations.  Since 
these  components  are  genera!  small,  this  coupling  is  expected  to  be  weak.  It  is 
difficult  to  formulate  the  theory  of  such  interactions,  because  of  the 
complexity  of  the  radiation  modes  of  a  fiber  whose  index  difference  is  not 
small.  Such  a  calculation  may  be  incorrect  in  principle  anyway,  because  the 
primary  coupling  to  radiation  modes  may  not  be  from  the  fundamental  mode,  but 
rather  through  a  "ladder"  excitation,  where  light  from  the  fundamental  mode 
scatters  into  higher  order  guided  modes,  which  in  turn  couple  more  strongly  to 
the  radiation  modes  than  does  the  fundamental  mode.  Thus  it  is  likely  the 
dominant  mechanism  when  direct  coupling  is  strongly  suppressed.  The  analysis 
in  this  case  is  complex,  so  that  it  is  difficult  to  estimate  the  radiation  losses, 
but  a  significant  fraction  of  the  light  lost  to  higher  order  guided  modes  may 
well  leak  out  of  the  fiber  One  consequence  of  such  a  higher  order  loss 
mechanism  would  be  a  length-dependent  radiation  loss  coefficient,  as  the  loss 
would  depend  on  the  mix  of  guided  modes  present  in  the  fiber.  A  length-resolved 
loss  measurement  would  be  necessary  to  determine  the  dominant  loss 
mechanism,  but  this  has  not  been  done  as  yet. 

Scattering  loss  measurements  have  been  made  by  a  method  described  in 
detail  in  Ref.  4.  The  results  obtained  at  Stanford  for  two  unclad  170)im 
diameter  sapphire  fibers  are  shown  in  the  Fig.  4.3.  One  of  the  fibers  was  grown 
without  feedback  control  of  the  diameter,  while  the  other  was  grown  with 
closed  loop  control.  Their  RMS  diameter  variations  are  2%  and  0.2% 
respectively  The  losses  are  quite  small,  on  the  order  of  10"^dB/cm,  and  are 
approximately  inversely  proportional  to  the  wavelength.  The  lowest  loss, 
measured  at  3.39pm,  is  equivalent  to  80  dB/km,  comparable  to  the  losses 
obrtained  in  early  glass  fibers.  The  loss  coefficient  of  the  diameter  controlled 
fiber  is  the  same  or  slightly  larger  than  that  of  the  open- loop  fiber.  This  can  be 
explained  by  the  characteristic  spatial  frequencies  of  the  diameter  variations  in 
the  closed-loop  fiber  which  are  higher  than  those  in  the  open- loop  fibers  and 
which  would  enhance  the  coupling  to  radiation  modes  in  the  closed-loop  fibers. 
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c.  Absorption  Losses 


The  absorption  losses  in  single  crystal  fibers  should  be  similar  to  those  int 
he  starting  rod  from  which  the  fiber  is  grown.  For  most  materials,  the  problem 
of  producing  fibers  with  low  absorption  losses  thus  reduces  to  the  problem  of 
finding  sufficiently  pure  starting  material.  Exceptions  to  this  rule  include 
materials  containing  a  volatile  component,  in  which  case  the  fiber  could  have  a 
different  composition  than  the  source,  and  materials  that  partially  decompose 
on  melting.  For  most  materials  of  interest  in  device  applications,  the 
absorption  is  small  compared  to  the  losses  due  to  guided  mode  coupling. 

The  measurement  of  absorption  losses  by  a  direct  insertion  loss 
measurement  is  not  appropriate  for  materials  with  an  absorption  that  is  less 
than  0.00 1cm'1,  because  of  the  difficulty  in  separating  the  small  reduction  in 
throughput  due  to  absorption  from  the  much  larger  reduction  due  to  reflection. 
The  simplest  technique  is  vacuum  calorimetry.  In  this  method,  the  sample  is 
illuminated  with  a  laser  beam  at  the  wavelength  of  interest,  and  the 
temperature  rise  is  measured  as  a  function  of  time.  The  temperature  rise  can 
be  related  through  simple  heat  transfer  calculations  to  the  amount  of  power 
deposited  in  the  fiber,  from  which  the  absorption  coefficient  can  be  extracted  if 
the  incident  laser  power  is  known. 

For  a  cylindrical  rod,  the  relation  between  the  absorption  coefficient  and  the 
temperature  rise  is; 

a=  nDhAT/P  42 

where  a  is  the  absorption  coefficient,  D  is  the  diameter  of  the  fiber,  H  is  the 
coefficient  of  heat  loss  from  the  surface  of  the  fiber,  aT  is  the  temperature 
rise,  and  P  is  the  incident  laser  power.  H  may  be  calculated  theoretically,  or 
extracted  from  the  exponential  decay  of  the  temperature  along  the  length  of  a 
fiber  heated  at  one  end.  Using  the  latter  method,  it  is  found  that  for  AI2O3  rods 

in  air,  h~  1  mW/Kcm2,  while  in  vacuum  h*0.3mW/Kcm2.  For  a  500  mm  diameter 
rod  we  have  that. 

a(cm_1)  =  0.16AT(°K)/P(mW)  43 

Stanford  has  measured  the  absorption  coefficient  in  "craquelle"  sapphire 
soui ce  material  42  At  51 5nm,  the  absorption  is  IxIO'^cm"'.  The  losses 
measured  in  fibers  grown  from  this  material  are  similar  to  those  for  the  source, 
as  would  be  anticipated  from  the  previous  discussion.  We  see  that  the 
absorption  and  radiation  loss  coefficients  (Fig.  43)  are  comparable,  so  both 
effects  must  be  considered 

LaserGenics,  in  a  program  supported  by  the  Army  Matei  lals  Technology 
Laboratory,  has  developed  an  extremely  sensitive  vacuum  calorimetry  technique 
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which  is  applicable  to  any  optical  material  as  well  as  optical  fibers.  We  have 
applied  the  technique  to  over  30  samples  of  LiNbOy  which  have  been  grown  by 

Crystal  i  prhnoiogy  inc  and  have  demonstrated  a  sensitivity  t nr  the 
meac.iirem^nt  of  thr*  absorption  coefficient  of  better  than  0  0003 

0  pA^nlty 


1  be  objective  of  tms  single-crystal  fiber  research  is  the  demonstration  of 
wffin^nr  nor* linear  optmai  devices  as  is  established  in  the  pre^iou^  sections, 
nonlinear  optical  devices  require  fibers  meeting  stringent  material,  structural 
and  processing  criteria  The  mode  confinement  of  guided  wave  structures  is 
useful  m  improving  the  mode  overlap  with  poor  quality  pump  sources,  for 
reducing  the  threshold  in  low-gain  materials  for  such  nonlinear  devices  as 
optical  parametric  oscillators,,  and  in  producing  very  high  gains  for  optica' 
conversion. 

The  scatter  losses  measured  in  fiber  lasers  have  been  shown  to  bo  in 


qualitative  agreement  with  the  predictions  of  Section  2.  Therefore,  we  can  have 
confidence  in  the  use  of  these  results  to  predict  losses  in  nonlinear  devices  in 
particular,  it  has  been  established  that  it  is  possible  to  propagate  a  single  mode 
in  highly  multimoded  structures,  even  without  the  benefit  of  a  low  index  of 
refraction  claodmo.  An  eyamp!e  of  the  output  mode  of  an  eadv  unclad 
monontnic  noyaG  finer  iaser^O  is  shown  in  Fig  4.4,  We  see  that  the  profile  is 

*  i 

close  tr>  Gaussian  More  recent  resultsM'  have  shown  more  than  at  tne 
output  power  in  the  fundamental  mode,  and  losses  on  the  order  of  0.01cm" '. 


ibis  level  of  performance  is  adequate  for  some  interesting  nonlinear  device? 
i  he  wort  on  i  iNnOy  fibers  has  focussed  on  tne  growth,  poling,  dadoing,  and 

characterization  studies.  Nonlinear  interactions  in  LiNbCh  fibers  have  been 


limited  fc  temperature  pnasematched  second  harmonic  generation  of  CM2  nrn 
ragtatior.  in  unclad  fibers,  with  efficiencies  comparable  to  these  of  bulk 
interactions  in  a  crystal  of  the  same  length. 

Early  studies  of  $H0  m  MgOLiNbOj  fibers*^6  suggested  that  lithium  was  lost 


during  growth,  presumably  through  volatilization  of  L12O.  Carrying  out  the 
growth  ib  a  wet  oxygen  atmosphere  substantially  reduced  this  problem,  but 
interest  in  the  precise  control  of  the  u/Nb  ratio  led  to  investigations  of  a 
wapor -phase  technique  for  Li  transfer,  first  suggested  by  Holman*46  as  a  means 
of  investigating  the  phase  diagram  of  L1NP03  The  method  involves  heating  the 


fiber  m  the  presence  of  a  powder  made  of  a  Li-rich  mixture  of  L i 3 f4bO 4/ L i N bO 3 


The  1 1  vapor  m  equilibrium  with  this  mixture  is  absorbed  by  the  lithium-poor 
fiber  until  tne  composition  of  the  fiber  reaches  the  phase  boundary  between  the 
1  inn '3  and  1  inn  13  =  Li3Nt'0«4  This  composition  is  essentially  stoichiometric. 


i.e .  the  U/Nb  ratio  is  close  to  one.  The  position  of  this  phase  boundary  is 
almost  independent  of  temperature,  so  the  choice  of  temperature  affects  only 


trie  rate  of  the  process,  not  the  end  result.  For  Li-poor  powders,  the  phase 
boundary  is  temperature  dependent,  and  the  composition  of  the  crystal  can  be 
controlled  by  the  choice  of  temperature. 

This  process  has  two  features  of  interest.  The  composition  of  the  crystal  is 
determined  by  thermodynamic  equilibrium,  and  becomes,  in  principle, 
independent  of  its  original  composition.  Thus,  it  is  possible  to  produce  crystals 
of  very  uniform  composition,  at  Li/Nb  ratios  other  than  the  congruently  melting 
one  Uniform  composition  is  important  for  efficient  nonlinear  interactions. 

The  other  important  feature  is  that  the  birefringence  of  lithium  niobate 
increases  with  the  Li/Nb  ratio.  Thus,  at  a  given  wavelength,  the  phasematchinq 
temperature  increases  with  Li/Nb  ratio,  so  that  one  can,  for  example,  produce  a 
composition  that  phasematches  above  the  annealing  temperature  for 
Photorefractive  damage.  Alternatively,  at  a  fixed  temperature,  the 
phasematching  wavelength  decreases,  so  that  it  may  be  possible  to  phasematch 
for  SH6  into  the  blue  at  room  temperature. 

Although  one  can  have  some  degree  of  control  of  the  quality  and  composition 
of  the  fiber  after  the  growth  of  the  fiber,  it  is  certainly  prefered  to  be  able  to 
grov;  fibers  of  high  optical  quality.  From  a  commercial  standpoint,  the  fewer 
steps  needed  to  produce  a  device  the  better  from  the  stand  point  of  yield  and 
therefore  cost.  During  the  course  of  the  program,  we  made  18  growth  runs  in 
our  facility  of  MgO.LiNbOj.  The  early  growth  runs  were  similar  to  those  grown 

at  Stanford  University  The  fibers  exhibited  a  significant  degree  of  browning 
and  nao  a  nigh  qpn^ity  of  defects  As  experience  was  gained  in  the  growth 
process,  rqp  quality  of  the  fibers  experienced  a  marked  improvement  in  optica) 
quality,  towards  tne  ^nd  or  the  program  we  were  able  t"  grow  MgO  L iNbOj 


fibers  of  very  good  optical  quality.  Some  examples  of  these  nb?rs  are  shown  m 
Fig.  4.5  Although  the  optical  quality  of  these  fibers  was  high  the  mechanical 
characteristics  'were  still  not  ideal  in  that  the  fibers  were  still  too  brittle  and 
therefore  too  fragile.  Annealing  the  fibers  markedly  improved  the  mechanical 
strength  of  the  fibers  but  we  were  not  able  to  eliminate  this  process  step 
Although  we  were  not  able  to  eliminate  the  annealing  process,  this  process 
for  high  quality  optical  fibers  is  rather  simple.  Our  biggest  disappointment 
came  in  the  area  of  poling  the  fibers  to  obtain  single  domain  fibers.  We  were 
not  able  to  maintain  the  necessary  stability  with  the  asymmetric,  heating 
technique  at  the  temperatures  required  to  achieve  single  domain  fibers.  We 
frequently  overheated  the  fiber  so  that  we  melted  through  it.  Combining  the 
electric  poling  w'lth  the  asymmetric  heating  technique  only  complicated  the 
process  The  only  process  that  seemed  to  appear  promising  was  the  combination 
of  uniform  heating  in  a  reaction  tube  and  electric  field  poling.  We  had  much 
better  control  in  this  process.  The  problem  with  our  experiments  with  this 
techmaue  was  that  at  this  time  the  furnace  we  had  was  only  capable  of  reaching 
1  iOO'C,  not  hiqh  enouqh  to  qet  close  to  the  melting  point  of  MqO.LiNbCh  at 


i  .v„vi  sinrp  topsp  pyperiments  were  dc»ne  w3  hav^  purchased  a  new  furnarp 


capable  of  reaching  1500°C.  We  intend  to  pursue  this  technique  in  the  future, 
however  it  must  be  noted  that  careful  control  of  the  temperature  is  required  as 
any  asymmetry  in  the  temperature  will  lead  to  the  distructlon  of  the  fiber.  Care 
must  also  be  used  in  holding  the  fiber  so  that  the  fiber  Is  in  close  contact  with 
the  electrodes.  Our  approach  to  this  problem  was  to  cut  a  fine  slot  in  a  high 
temperature  dielectric  such  as  sapphire  which  holds  the  two  high  temperature 
wire  electrodes  and  the  MgOliNbCh  fiber.  We  are  confident  that  this  technique 

wiil  succeed  after  further  development 

e  Conclusions 

wp  have  found  that  MgOUNbCw  can  be  grown  In  fiber  form  with  very  good 

optical  quality.  Although  the  mechanical  strength  is  not  what  is  required  for 
the  as-grown  fibers,  the  fibers  can  be  strengthened  with  a  rather  simple 
annealing  process  The  greatest  difficulty  was  encountered  in  attempting  to 
pole  the  fibers  in  order  to  obtain  single  domain  a-axis  fibers.  However,  a 
process  utilizing  both  heating  and  electric  field  poling  does  have  promise  for 
achieving  the  single  domain  fibers  required  for  nonlinear  devices. 
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RELATIVE  INTENSITY  (%) 


Figure  4.4.  Mode  profile  of  an  early  unclad  Nd'.YAG  fiber  laser 
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Figure  4.5.  Fibers  of  MgO:LiNbOs  grown  by  the  LIIPG  technique. 
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